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Abstract 

The muon energy spectra of the quasi-elastic and 1-pion production events in a LBL 
experiment, hke K2K, are predicted to follow closely the neutrino energy spectrum, 
with downward shifts of the energy scale by < Q2 > /2M and (< > +M\ - 
M^)/2M respectively. These predictions seem to agree with the observed muon 
spectra in the K2K nearby detector. The corresponding muon spectra in the far- 
away (SK) detector are predicted to show characteristic spectral distortions induced 
by oscillation. Comparison of the predicted spectral distortions with the observed 
muon spectra of the 1-Ring and 2-Ring muon events in the SK detector will help to 
determine the oscillation parameters. The results will be applicable to other LBL 
experiments as well. 
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Recently the KEK to Kamioka long-baseline neutrino oscillation experiment (K2K) 
has published its result [1], which confirms the existence of i/^ oscillation as seen in the 
Supcr-Kamiokande (SK) atmospheric neutrino data [2]. It uses an accelerator-produced 
i^fj^ beam, which has a peak energy of 1.2 GcV and a spread of about 1 GeV on either side - 
the full width at half maximum being ~ 1 GeV [1, 3]. The beam passes through a nearby 
neutrino detector (ND), comprising of a 1 kiloton water Cerenkov (IKT) and a fine- 
grained detector (FGD) system, and then travels a distance of 250 km to enter the SK 50 
kiloton water Cerenkov detector. The charged- and neutral-current event rates measured 
in the ND are used to predict the corresponding event rates at the SK detector in terms 
of the oscillation parameters. They report 56 fully contained (FC) events against the 
expectation of 80 without oscillation. The observed deficit agrees well with the neutrino 
mass and mixing angles deduced from the atmospheric neutrino oscillation data [2] 

sin^ 29 ~ 1 and Am^ ~ 3 x lO'^eV^ . (1) 

We estimate that the QE, 1-pion production and multi-pion events constitute 35%, 50% 
and 15% at E,, — 1.2 GeV. In this region the A(1232) resonance contribution ac- 
counts for 80 — 90% of the 1-pion events [4], with the remainder produced from the 
Pii(1440), 5ii(1535) and L>i3(1520) resonances. 

In this note we use kinematic considerations to predict the muon energy spectra of 
the QE and 1-Pion resonance production events which constitute the bulk of the charged- 
current z/^ scattering events in the K2K experiment. These predictions can be checked 
with the observed muon energy spectra from the ND. We also present the distortion of 
these muon spectra due to oscillation, which one expects to see at the SK detector. 
Comparison of the predicted muon spectra with those of the observed QE and 1-Pion 
events at the SK detector will be very useful in determining the oscillation parameters. 

Since the A resonance width of 120 MeV is very small compared to the spread of the u 
beam energy, one can make a narrow resonance approximation for A. A common property 
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of QE and narrow resonance production processes is a highly restricted kinematics, i.e., 

W^^M^ -Q^ + 2ME - 2ME' , (2) 

where W represents the nucleon mass M for QE and Mr for resonance production. The 
other variables are as usual the neutrino energy {E), the muon energy {E') and the 4- 
momentum squared {Q^) transferred between them. We can approximate the production 
cross-section for the narrow resonance as 

-j^ = |,,(^^-Ml). (3) 

where we have replaced the Breit-Wigner factor by the 5-function. Note that this relation 
holds exactly for QE scattering with Mji = M. Thus we get from eqs. (2) and (3) 

_^^ = ^S(E'-E+^M^) (4) 

where AM^ = — M^, which is for QE scattering. Convoluting this cross-section 
with the neutrino beam spectrum f{E) gives the flux averaged cross-section 

,) . . E,) HE ^E,).E,^E'+^. (5) 

Moreover to a first approximation one can integrate over to write 

/ f{E) ^ d£; = a^{E' + ^^^q^) m + ^M^^) , (6) 

where < >~ 0.3 GeV^ for both QE and A production and the neutrino energy E — 
— 3.5 GeV for the K2K experiment [1]. It assumes the effective range of Q'^ integration 
to be small, so that over the corresponding neutrino energy range Ai? = A(5^/2M the 
f{E) and a^[E) can be approximated by their mean values. Thus we see from eq. (5) or 
(6) that for both QE and narrow resonance production the shape of the muon spectrum 
is predicted to follow the shape of the incident neutrino spectrum, with a downward shift 
of the energy scale by ^^2m^^ ■ notes an interesting analogy between this result and 
the quark parton model, where the structure function in deep inelastic scattering follows 
the shape of the parton energy distribution within the nucleon. 
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Of course the shape of the muon spectrum is not exactly identical to the shape of 
the neutrino spectrum, because it is modulated by the QE (Resonance production) cross- 
section. The QE (A production) cross-section shows a threshold rise over E = — 1 GcV 
(0.4 — 1.4 GeV) and becomes flat thereafter. This will be reflected in steepening of the 
muon spectra at the low energy end, as we shall see below. On the other hand the shape of 
the muon spectra on the higher energy side of the peak will be practically identical to that 
of the incident neutrino spectrum for both QE and A production processes. Therefore the 
relative size of the SK to the ND cross-sections over this region provides a direct measure 
of the spectral distortion and hence the underlying oscillation parameters. 

These are simple but very useful results, particularly since A accounts for about 
80 — 90% of the 1-Pion production cross-section and the QE plus 1-Pion production 
together account for 85% of the charged-current i/^ cross-section at the K2K experiment. 
It means that one can predict the muon energy spectra of the QE and 1-Pion events 
in terms of the neutrino energy spectrum from general kinematic considerations without 
practically any model dependent ansatz. These predictions can be verified with the muon 
spectra of the QE and 1-Pion events observed by their ND. Even more importantly one 
can then predict the distortions in the corresponding muon energy spectra at the SK 
detector, induced by oscillation. Thus the predictions corresponding to eqs. (5) and 
(6) for the SK detector are 

^^^^^^ WdE^ ^ dQ^^^^^ ^^^^^ " > (7) 

with = + and 

fsK{E) -^dE^ a^{En) f{En) [1 - sin^ 29 sm'C-^^)] , (8) 

with Eji — E' + ^^'^ ^ . Comparing these predictions with the observed muon energy 
spectra of the QE and 1-Pion events of the SK detector will test the spectral distortion 
due to oscillation and determine the oscillation parameters. 
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Let us first present a simple analytic form of QE and A production cross-sections for 
the ND and SK detector corresponding to eqs. (6) and (8) respectively. For this purpose 
we have approximated the neutrino energy spectrum by the Lorentzian 

f{E) = — J ^ , Eo = 1.2 GeV , 7 = 0.6 GeV , (9) 

which agrees with the shape of the K2K neutrino energy spectrum quite well. The nor- 
malization factor of 1.25 ensures that the total flux is normalized to 1 when integrated 
over the experimental range of £^ = — 3.5 GeV instead of — oo to oo. Figure 1 shows 
a comparison of the K2K neutrino energy spectrum with the parametrisation of eq. (9). 
Substituting this in eqs. (6) and (8) gives 

//(.) g , (10) 



and 



1 - sm 20 sm | _ AMH<g^> 



2M 
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Finally we have approximated the QE and A production cross sections by 

a^{E) = En X 4.5 ib , En = uim{E, 1) GeV , (12) 
a^(E) = (Ed - 0.4 GeV) x 5.5 fb , Ed ^ mm{E, 1.4) GeV . (13) 

They correspond to linear threshold rise for the QE (A production) cross-section between 
0-1 (0.4-1.4) GeV, which describes the cross-sections quite well. The normalizations 
correspond to the average cross-section per nucleon for a H2O target [5]. 

For the sake of completeness we have also calculated the Pii(1440) and 5'ii(1535) 
contributions using the narrow resonance approximation. Of course it may not work 
so well for the Pii(1440), which has a width of 350 MeV. But since these resonance 
contributions to the 1-Pion production cross-section are small, the resulting error will be 



insignificant. Although there are significant differences between the parametrisations of 
ref. [5, 6] and ref. [7] for the individual resonances, the cumulative contributions are very 
similar for the two sets of parameters. We have used the parametrisation of ref. [5, 6] for 
these resonance production cross-sections. 

Figure 2a shows the predicted muon energy spectra for the QE and 1-Pion events, the 
latter being dominated by the A contribution (shown separately). One can clearly see 
the downward shift of the peak from 1.2 GeV hy < > /2M ~ 0.15 GeV for the QE 
and < >)/2M ~ 0.4 - 0.5 GeV for the A contribution. One can also see the 

reflection of the threshold rise from eqs. (12) and (13), particularly for the A contribution. 
On the other hand the muon energy spectra to the right of the peak closely follow the 
shape of the neutrino energy spectrum of Fig. 1. Both the curves agree reasonably well 
with the more exact calculation of QE and 1-Pion production cross-sections, except for a 
small range near E' ~ 0.4 GeV as discussed below. One can see them to be in general 
agreement with the observed muon spectra of QE and non-QE events of the K2K ND 
data [1]. But a quantitative comparison will need to include things like the detection 
efficiency and experimental cuts, which is beyond the scope of the present work. 

Figure 2b shows the corresponding muon energy spectra of the QE and 1-Pion events 
for the far away detector, predicted by eqs. (11, 12, 13) along with a small contribution 
from the higher resonances. One can clearly see the distortion of the muon energy spectra 
due to oscillation, where we have used the oscillation parameters of eq. (1). They 
are again in reasonable agreement with the more exact results presented below. It will 
evidently be very useful to compare these predictions with the observed muon energy 
spectra of QE and 1-Pion events from the SK detector which is not available to us. Again 
a quantitative comparison will require taking into account the detection efficiency of the 
muon and the pion, which can best be done by the K2K Collaboration. 

We shall now present the results of an exact calculation. It is based on numerical 
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integration of eqs. (5) and (7) for QE scattering, while for the resonance production it 
uses the corresponding formulae including the resonance widths. Apart from the leading 
contribution from the A(1232) resonance, it includes subleading contributions from the 
Pii(1440) and 5*11(1535) resonances. Wc estimate the contribution from the still higher 
resonances along with the nonresonant background to be no more than 5 — 10% of the 
1-Pion production cross-section at this energy. Therefore the accuracy of our prediction 
should be as good as that of the K2K experiment. The parametrisation of the QE and 
the resonance production cross-sections used in this analysis can be found in [5, 6]. 

For the excitation of the A- resonance we use the formahsm of [8, 9] with two distinct 
form factors for the vector and axial currents. The vector and the axial vector form factors 
are modified dipoles [10]. For the Pu and resonances we use the form factors from 
[11]. 

Fig. 3a shows the predicted muon spectra for the QE and 1-Pion production processes. 
One can clearly see the downward shift of the peak from 1.2 GeV by < > /2M for the 
QE and (AM^-F < >)/2M for the A contribution, dominating the 1-Pion production 
process. The steepness of the muon spectra at the low energy end reflects the threshold 
rise. On the other hand they closely follow the shape of the neutrino energy spectrum on 
the right side of the peak. All these features are in general agreement with the approximate 
results shown in Fig. 2. One sees however a clear broadening of the muon spectra below 
1 GeV compared to the latter. In particular the exact calculation for QE shows a clear 
excess over the approximate one near E' ~ 0.4 GeV as mentioned earlier. This is because 
the assumption of a small range of integration made in deriving cq. (6) from eq. (5) does 
not hold here, as the rapid increase of the flux factor upto E = E' + Q'^/2M = 1.2 GeV 
stretches the effective range of Q'^ upto 1.5 GeV^. Note that below E' = 0.4 GeV the 
range is restricted by the kinematic constraint Q'^ < AME'"^ / (M — 2E'). So the rapid fall 
of the spectrum with E' simply reflects the fall of the flux and the threshold factors with 
E. 
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The predicted muon spectra of Fig. 3a agree reasonably well with the corresponding 
spectra of the K2K ND [1] for both the QE and the non-QE parts. In particular one 
can compare the predicted QE spectrum with their simulated QE spectrum shown in Fig. 
1 of ref. [1]. Their Fig. la and c show separately the QE muon spectra for the 1-Ring 
muon (IR//) sample of the IKT and the QE enhanced sample of the FGD respectively. 
The two play complementary roles in covering the complete muon energy range, as the 
IKT and the FGD have high efficiencies at < 1 GeV and E' > 1 GeV, respectively 
[1]. One can not compare our predicted muon spectra with these figures quantitatively 
without folding in these efficiency factors, which are not available to us. But there is good 
qualitative agreement between the predicted QE spectrum of our Fig. 3a with their Fig. 
Ic at £" > 1 GeV and Fig. la at £" < 1 GeV. While the former shows the position of 
the peak and the shape of the spectrum to the right, the latter shows the broadening of 
spectrum down to E' ~ 0.4 GeV as discussed above. Similarly one sees good agreement 
between the predicted muon spectrum of our Fig. 3a for 1-Pion events with the non-QE 
spectra of their Fig. lc,d at £" > 1 GeV and Fig. la at £" < 1 GeV. Thus one has a 
simple and robust prediction for the shape of the muon spectrum in terms of the neutrino 
spectrum not only for the QE events but also for the 1-Pion events, which dominate the 
inelastic events. 

Fig. 3b shows the corresponding muon energy spectra of the QE and 1-Pion events 
for the SK detector, predicted by eqs. (1) and (7). Again the curves are in reasonable 
agreement with the approximate results of Fig. 2b. One can clearly see the distortion 
of the muon energy spectrum due the oscillation. They should be compared with the 
observed muon energy spectra of the 1-Ring and 2-Ring muon events at the SK detector, 
after taking into account the pion detection efficiency. We hope such a comparison will 
be done by the K2K collaboration. 

To estimate the pion detection efficiency one needs to take into account the nuclear 
absorption, Pauli blocking and charge exchange taking place during the rescattering of the 



7 



produced pions. We have included these effects following the prescription of ref. [12] . Since 
the dominant contribution to 1-Pion production process comes from resonance production 

on Oxygen, we have evaluated the effects of nuclear absorption and rescattering on the 
produced pion for this case. The relevant subprocesses are up //."pTr"*", un fi^nir^ 
and un ^"jm'^ with relative cross-sections 9 : 1 : 2 for the dominant contribution from 
A. 

Fig. 4 shows the effects of nuclear corrections on the produced 7r+ and 7r° spectra from 
these processes for the nearby detector averaged over the neutrino spectrum. The results 
are very similar for the SK detector. Nuclear rescattering effects result in enhancing the 
tt" events at the cost of the dominant tt^ component. But taken together we see a nearly 
20% drop in the rate of 1-Pion events due to nuclear absorption of the produced pion. 
Moreover about 10% of the remaining events corresponds to the pion momentum being 
less than the Cerenkov threshold of 100 MeV. Therefore one expects about 70% of the A 
events to give a detectable pion ring at the SK detector while the remaining 30% appears 
hke a QE event. Adding the latter to the 35% of genuine QE events would imply that 
about 50% of the CC events will appear QE-likc at the SK detector. This seems to agree 
with the observation of 29 1-Ring muon (IR//) events at the SK detector out of their 56 
FC events. It also agrees with their Monte Carlo simulation of 44 IRji events out of the 
total of 80 FC events in case of no oscillation. Thus the observed muon energy spectrum 
of the 29 IRji events should be compared with a weighted average of the predicted spectra 
for QE and 1-Pion events. Alternatively the observed muon energy spectrum of the sum 
of 1-Ring and 2-Ring muon events could be compared with the predicted spectrum of 
the sum of QE and 1-Pion events. Although a part of the 2-Ring events may come from 
multi-pion production, the resulting error may be small since multi-pion events at the ND 
constitute only ~ 15% of CC events. 

It should be mentioned here that the K2K Collaboration has reconstructed the neu- 
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trino energy for the 29 IR// events using the kinematic relation 

ME' - ml/2 

E = ^ 14 

M -E' + P^ cos 0^ ^ ^ 

and used it for testing the spectral distortion due to neutrino oscillation [1]. While this 
relation holds for genuine QE events, it is not possible to reconstruct the neutrino energy 
for an inelastic event using only the muon parameters. As we have seen above, about 
l/3rd of the 29 IR/x events may come from 1-Pion production. Hence we can see no 
reliable way of reconstructing the neutrino energy for this IR/x sample on an event by 
event basis. On the other hand the muon energy is a directly measurable quantity for each 
event. Therefore it seems to us to be a better variable for testing the spectral distortion 
phenomenon compared to the reconstructed neutrino energy. 

Let us note finally that the above formahsm is apphcable not only to K2K but also to 
future long basehne experiments like J2K, MINOS and the CERN-Gran Sasso experiments 
which plan to use low energy beams [3, 13]. Therefore it will be very useful to extend 
this analysis for the beam energy spectra and the target nuclei of these experiments. 
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Figure 1: K2K neutrino energy spectrum. The solid line is the exact spectrum normalised to 
unit area. The dashed line shows the approximated spectrum by the Lorentzian in Eq. (9). 
We also show for comparison as dotted line the same Lorentzian, but with normalisation 1.4 
instead of 1.25. 
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Figure 2: Approximate prediction of the muon energy spectra for the (a) Nearby and (b) SK 
detectors of the K2K experiment. The Quasi-Elastic and the 1-pion production cross-sections 
are shown along with the A contribution to the latter. 
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Figure 3: Exact predictions of the muon energy spectra for the (a) Nearby and (b) SK detectors 
of the K2K experiment. The Quasi-Elastic and the 1-Pion production cross-sections are shown 
along with the A contribution to the latter. 
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Figure 4: The momentum distribution of the decay pion for resonance production by charged 
current in Oxygen before (dotted) and after nuclear correction (solid). The dashed line shows 
only the effect of Pauli blocking. 
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